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The effect of out-of-plane orientational freedom on the orientational ordering properties of a 
monolayer of hard ellipsoids is studied using Parsons-Lee scaling approach and replica exchange 
Monte Carlo computer simulation. Prolate and oblate ellipsoids exhibit very different ordering 
properties, namely, the axes of revolution of prolate particles tend to lean out, while those of oblate 
ones prefer to lean into the confining plane. The driving mechanism of this is that the particles try to 
maximize the available free area on the confining surface, which can be achieved by minimizing the 
cross section areas of the particles with the plane. In the lack out-of-plane orientational freedom the 
monolayer of prolate particles is identical to a two-dimensional hard ellipse system, which undergoes 
an isotropic-nematic ordering transition with increasing density. With gradually switching on the 
out-of-plane orientational freedom the prolate particles lean out from the confining plane and a 
destabilisation of the in-plane isotropic-nematic phase transition is observed. The system of oblate 
particles behaves oppositely to that of prolates. It corresponds to a two-dimensional system of hard 
disks in the lack of out-of-plane freedom, while it behaves similar to that of hard ellipses in the 
freely rotating case. 


I. INTRODUCTION 


Phase behavior of non-spherical hard bodies with their 
centers of mass confined in planar geometry is receiving 
considerable attention due to the recent development of 
the preparation of colloidal particles with various shapes 
and new experimental techniques. Nowadays it is pos¬ 
sible to prepare colloids with several geometrical shapes 
such as cubes, polyhedrons, octopods, ellipsoids and he¬ 
lices [1-6]. The anisotropic colloids can be confined at 
the interfaces [7-11], between two parallel solid walls 


* Electronic address: vargasz@almos.uni-pannon.hu 
1 Electronic address: yuri@math.uc3m.es 
^Electronic address: enrique.velasco@uam.es 
§ Electronic address: gustavo.bautista@uacm.edu.mx 
^Electronic address: godriozo@azc.uam.mx 


[12,13], at the bottom of the sample holder [14], at a 
substrate surface [15] and into a lamellar matrix of sur¬ 
factants [16,17]. The confinement can be so strong that 
even colloidal monolayers can be realized experimentally. 
Ordering properties of two-dimensional and quasi two- 
dimensional (q2D) non- spherical colloids has been the 
subject of several experimental and theoretical studies 
[17-36]. The reason for this is that the nature of two- 
dimensional (2D) nematic ordering is quite different from 
the three-dimensional (3D) one. It shows only quasi-long- 
range orientational order with algebraically decaying ori¬ 
entational correlations and the ordering transitions be¬ 
tween isotropic and nematic phases are first order or con¬ 
tinuous through a Kosterlitz-Thouless disclination un¬ 
binding type mechanism [18, 37]. Strictly 2D colloidal 
systems cannot be examined experimentally, because the 
out-of- plane orientational and positional freedoms are 
always present to some extent. Therefore it is desir¬ 
able to extend the theoretical studies in such directions. 
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where the extra orientational and positional freedoms are 
present. Along this line, the ordering properties of mi¬ 
crotubules confined in a thin slit have been modeled as 
hard spherocylinders placed between two planar walls in 
[28]. In agreement with the experiment it has been found 
that the isotropic-nematic transition density increases 
with the wall separation [28]. Ellipsoidal shaped colloidal 
particles are gaining widespread applications due to de¬ 
velopment of the stretching techniques in preparation of 
monodisperse prolate and oblate ellipsoids from spherical 
latex particles [4,38]. They can be also confined into pla¬ 
nar geometry and create a monolayer to study ordering 
and glassy behavior of the q2D ellipsoid systems [12,13]. 

In this paper we examine the orientational ordering 
properties of q2D hard ellipsoid systems, where parti¬ 
cles are allowed to rotate out the confinement plane to 
some extent, while it is assumed in first approximation 
that the centers of the particles are always in the same 
plane. Switching on the out-of-plane orientational free¬ 
dom by gradual increase (decrease) of the limiting polar 
angle (Oc) for oblate (prolate) shapes, it is possible to 
make a link between strictly 2D hard ellipse (hard disk) 
systems for prolate (oblate) shapes and q2D freely ro¬ 
tating prolate (oblate) ellipsoids systems. We show that 
the additional out-of-plane orientational freedom changes 
substantially the orientational ordering and the transi¬ 
tion properties of both oblate and prolate shaped ellip¬ 
soids. To maximize the available free area on the confin¬ 
ing surface, the axis of revolution of the prolate particle 
leans out from the plane, while the oblate particle leans 
into the plane. As a consequence, the freely rotating 
prolate ellipsoids resemble hard disks at high densities, 
while those of oblate ellipsoids behave similarly to hard 
ellipses. The isotropic-nematic phase transition of hard 
ellipses corresponds to a planar nematic-biaxial nematic 
phase transition of the oblate ellipsoids. Here we note 
that our present work can be considered as an extension 
of our previous studies done for monolayers of uniaxial 
and biaxial hard particles using restricted-orientation ap¬ 
proximation [39,40]. 

The paper is organized as follows. The molecular 
model and the details of the confinement are presented 
in Sec. II. Sec. Ill is devoted to the Parsons-Lee the¬ 
ory of q2D hard ellipsoid fluids, where we show how to 
determine the equilibrium free energy, surface coverage, 
and order parameters of the uniaxial and biaxial nematic 
phases. Technical details of the replica exchange Monte 
Carlo simulation method are given in Sec. IV. The order 
parameters, the surface coverage, the equation of state 
and the phase diagram of the system hard ellipsoids are 
presented in Sec. V. Finally, some conclusions are drawn 
in Sec. VI. 


II. MOLECULAR MODEL 


In this work we examine the q2D system of hard el¬ 
lipsoids, where the particles are allowed to rotate freely 


Allowed orientations for prelates: d^<,d <.n-d^ 



Z 

y Allowed orientations for oblates: 



FIG. 1: Schematic representation of the system of confined 
hard ellipsoids: prolate ellipsoids (upper panel), oblate ellip¬ 
soids (lower panel). The particles are confined both position- 
ally and orient at ionally. The center of mass of the ellipsoid is 
always in the XY plane, which is denoted by the dashed line, 
while the ellipsoid can rotate freely in azimuthal angle (0) 
and is restricted by the polar angle {Oc)- The arrows indicate 
the direction of the main symmetry axis of the ellipsoids. 


in a restricted region of the solid angle, while the cen¬ 
ters of particles can move only on a two-dimensional XY 
plane. The shape of the particles can be both prolate 
and oblate, i.e. the aspect ratio or shape anisotropy ( 
K = cr||/cr_L, where and a± are the lengths along and 
perpendicular to the axis of revolution of the ellipsoid, 
respectively) can be either larger than one (/^ > 1 , pro¬ 
late shape) or between zero and one {0 < k < oblate 
shape). The orientational restriction takes place only in 
the polar angle (6>), which is measured from the axis of 
revolution of each particle to the normal of the confin¬ 
ing XY plane. While the particles can rotate freely in 
the azimuthal angle (0 < 0 < 27r), the allowed range of 
the polar angle is different for prolate and oblate ellip¬ 
soids, namely, the range of the polar angle is given by 
Oc < 0 < IT — Oc for prolate shapes, while 0 < 0 < Oc and 
7T — Oc < 0 < 7T intervals are allowed for oblate ones (see 
Fig. [^. We can tune Oc between 0 and 7r/2. In the case of 
prolate ellipsoids (/^ > 1) 6>c = 0 corresponds to the con¬ 
fined system of freely rotating particles, while Oc = 7r/2 
gives the two-dimensional system of hard ellipses. The 
situation is different for oblate ellipsoids (0 < /^ < 1), 
because the hard disk limit is given by = 0 ? while the 
freely rotating system of confined oblate ellipsoids corre¬ 
sponds to 6>c = 7r/2. The interaction between ellipsoids is 
purely hard, i.e. the pair potential between two particles 
is given by 


u(ri2,a>i2,wi,UJ2 



00 , 

0 , 


for ri 2 < ( 7 (u;i 2 ,a>i,u; 2 ), 

otherwise, ^ ^ 


where ri 2 is the center-to-center distance, 
u)i 2 = (cos 012 , sin 012 , 0 ) is the unit vector 
connecting the centers of the two ellipsoids, 
uji = (sin Oi cos (pi , sin Oi sin <pi , cos Oi) is the orienta¬ 
tional unit vector of particle i {i = 1 , 2 ) and cr(-• •) is 
the distance of closest approach. Note that the third 
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component of C 1 J 12 is always zero, which ensures that 
the centers of the particles are always in the XY plane. 
The distance of closest approach between two ellipsoids 
[cr(- • •)] is approximated in the manner proposed in our 
previous works [41,42]. 

III. PARSONS-LEE THEORY 

To describe theoretically the orientational ordering 
properties of the monolayer of hard ellipsoids we derive 
our working equations from the well-known Parsons-Lee 
(PL) theory of hard bodies [43,44], which proved very 
successful in determination of the equation of state and 
the transition properties of IN phase coexistence of non- 
spherical hard body fluids both in two [23,24] and three 
dimensions [45-47]. Here we present only the important 
ingredients of the theory, specific equations for the ellip¬ 
soid monolayer and some technical details. It is custom¬ 
ary to deal with the free energy of the system, which 
is the sum of ideal and excess free energy terms, i.e. 
F = Fid+Tex • The ideal term can be determined exactly 
from 

= logp-l + J dujf{u})logf{u}), (2) 


Y cri sin^ ^ where 0 is the polar angle 

of the ellipsoid. This means that the monolayer of the 
hard ellipsoids can be visualized as a multicomponent 
mixture of hard ellipses on the XY plane. Using the di¬ 
mensions of the hard ellipse, the intersected area of the 
hard ellipsoid with the XY plane can be obtained from 



This polar angle dependent area together with the orien¬ 
tational distribution function allows us to determine the 
surface coverage of the plane by ellipses, which is given by 
T] = p(a), where the average area of the intersected ellipse 
can be determined from (a) = J duja{0) Note that 

the surface coverage gives back the packing fraction of 
hard disks in the hard-sphere limit, because (j\\ = (j^ = D 
and the distribution function is a constant. In PL ap¬ 
proach the area of hard disks and the average area of 
hard ellipses are assumed to be the same (ukd = 
which implies also that 77 HD = V due to one-to-one corre¬ 
spondence between the two systems. Based on the con¬ 
cept of Parsons and Lee we can now calculate the excess 
free energy of the monolayer of hard ellipsoids using the 
thermodynamic properties of hard disks as follows 


where {3 = Ijk^T is the inverse temperature, N is the 
number of particles, p = N/A , A is the area of the 
plane, uj = {4>,0) is the collection azimuthal and polar 
angles, duj = d(f)dO sin 0 and /(cj) is the normalized ori¬ 
entational distribution function (f dujfiuj) = I). The 
ranges of the azimuthal and polar angles in the integra¬ 
tions have been already given in Sec. II. The excess free 
energy contribution can be obtained approximately with 
the mapping procedure from the actual system to a ref¬ 
erence one, where the second virial coefficient and the 
excess free energy can be obtained with good accuracy. 
In our case we choose the system of 2D hard disks as a 
reference system, because our confined ellipsoid system 
is q2D. The second virial coefficient of the hard disks is 
given by where D is the diameter of the 

hard disk. A luckily simple and accurate expression can 
be derived for the excess free energy of hard disks, too, 
using the scaled particle theory [48] 

= - In (I - t^hd) + , (3) 

1 - ^HD 

where 77 HD = is the packing fraction of the hard 

disk and ukd = /2 is the area of the hard disk. 

To perform the mapping procedure from the hard ellip¬ 
soids into hard disks we introduce the surface coverage 
(or packing fraction) of the hard ellipsoids, because the 
intersection of the ellipsoid with the XY plane is an el¬ 
lipse with characteristic lengths depending on the ori¬ 
entation of the ellipsoid. These characteristic lengths 
of the ellipse can be calculated easily using simple ge¬ 
ometry. One can derive that a±^e = and cr|| g = 



/3Fex /3Fex 

N ^ N SfiD’ w 

where the second virial coefficient of hard ellipsoids re¬ 
stricted to XY plane is given by 

-S™ = ^ J duJif{LOi) J dW2/(W2)^excl(Wl,W2). (6) 

This equation contains the excluded area between two 
ellipsoids, which can be obtained from the distance of 
closest approach (a) as follows 

I 7^^ 

Aexcl(wi,W2) = - y d<?!'i2Cr^(wi,W2,Wi2). (7) 

Using Eqs. (3)-(6) and t^hd = V condition the excess free 
energy of the system becomes 


N 


-ln{l-ri) + A— 

1 — T] 


J duJif{uJi) f da;2/(^2)^excl(^l,^2) 
4 f duja{0)f{uj) 


( 8 ) 


The sum of Eqs. (2) and (8) constitutes our density 
functional equation to determine the equilibrium orien¬ 
tational distribution function [f{oj)] and the total free 
energy density at a given density {p = N/A). The 
minimization of the free energy functional with respect 
to /(cj) gives us the Euler-Lagrange equation for the 
equilibrium f{uj). Note that the minimization must be 
carried out by maintaining the normalization condition 
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(f dujf{uj) = 1) . Although the minimization procedure 
is very simple and straightforward, we do not present 
the equation of /(cj) , since the resulting equation is too 
long. Once the equilibrium /(cj) is determined with a 
standard iterative method, the free energy density can 
be obtained by the substitution of the resulting /(cj) 
into the free energy functional (sum of Eq. (2) and (8)). 
We determine the pressure from the free energy using 
P = p^d{F /N) /dp and examine the orientational order¬ 
ing properties of the monolayer by the standard uniaxial 
and biaxial order parameters. These are defined as 

S = {P2) = J dujf{uj)P2{e), ( 9 ) 

and 

A = y* (ico’/(co’) sin^ 6>cos(20), (10) 

where P 2 {x) = 3x^/2 — 1/2 is the second order Legendre 
polynomial. The uniaxial order parameter {S) is posi¬ 
tive for out-of-plane ordering (0 < S' < 1), while it is 
negative for in- plane ordering (—1/2 < S < 0). The 
biaxial order parameter (A) can be very useful in finding 
in-plane orientational ordering transitions, because A is 
zero for in-plane complete disorder, while it is nonzero for 
in-plane order. S is practically a three-dimensional order 
parameter of bulk ellipsoids, while A is the correspond¬ 
ing two-dimensional one of bulk ellipses when Oc = 7r/2 
. Since our system is q2D, we need both of them. To 
see the effect of out-of-plane orientational freedom and to 
which extent the particles are tilted off from the confining 
plane, we determine the average aspect ratio of the effec¬ 
tive ellipse system, which is defined as /^eff = (<^||,e/<^-L,e)- 
In Sec. V we plot the density and the pressure in dimen¬ 
sionless units, which are defined as p* = Na‘j_/A and 
P^ =13 Pal. 

IV. REPLICA EXCHANGE MONTE CARLO 
SIMULATION 

As for the 2D hard ellipse study [35], we are implement¬ 
ing the replica exchange Monte Carlo technique [49-51]. 
This is done to avoid, as far as possible, the inherent 
hysteresis associated to transitions [52]. The method is 
based on the definition of an extended ensemble with par¬ 
tition function Qext = 11^=1 where Qi is the partition 
function of ensemble i. rir ensembles are considered, and 
Ur replicas are employed to sample the extended ensem¬ 
ble, each one at each ensemble. Defining Qext allows 
introducing swap trial moves between any two replicas, 
whenever the detailed balance condition is satisfied. In 
our case it is convenient to expand isobaric-isothermal 
ensembles in pressure [53]. This is so since we are study¬ 
ing hard particles. Hence, the partition function of the 
extended ensemble reads [53,54] 

Ur 

Qext = (11) 

i=l 


where QntPi is the partition function of the isobaric- 
isothermal ensemble of the system at pressure Pi and 
temperature T. N particles are taken into account at 
each ensemble. A standard implementation is used to 
sample the NTP i ensembles. This implies independent 
trial 2D displacements, 3D rotations of single ellipsoids, 
and area changes of the simulation cell. In case of hav¬ 
ing confining planes, 3D rotations of single ellipsoids are 
constrained by their presence. We are also accounting 
for non-orthogonal parallelogram cells and so, additional 
trial changes of the angles and relative length sides of 
the cell lattice vectors are included. The following accep¬ 
tance rule is set [53] 

-Prm,acc = min{l,exp [13 {Pi - Pj) {Ai - Aj)]} , (12) 

where Ai — Aj is the area difference between replicas 
i and j . Adjacent pressures must be close to provide 
swap acceptance rates over 0.1. Simulations are started 
from a packed triangular arrangement of spheres which 
are elongated in the direction normal to the plane by a 
factor K . In case of oblates, a stretching factor k. is also 
applied in a certain in-plane direction. Conversely to the 
stretching of spheres in a 3D cell, this procedure leads to 
the largest packed arrangement of spheroids in a plane 

[55] . A stationary state is reached faster by decompress¬ 
ing packed cells than by compressing lose random config¬ 
urations [52]. We perform the necessary trial moves to 
observe a stationary state. At this stage we adjust maxi¬ 
mum displacements to produce acceptance rates close to 
0.3. We also relocate all pressures, initially set by follow¬ 
ing a geometric progression with the replica index, to ob¬ 
tain similar swap acceptance rates for all pairs of adjacent 
ensembles. Next, we perform 4 x 10^^ sampling trials for 
fixed maximum particle displacements, maximum rota¬ 
tional displacements, maximum area changes, and maxi¬ 
mum changes of the lattice vectors. Verlet neighbor lists 

[56] are used to improve performance. We set N ^ 400 
ellipsoids and as a function of the pressure range to be 
covered. N ^ 400 is sufficiently large in view of Xu et. 
al. analysis of system size effects [33]. More details on 
the employed methods are given in previous works [52]. 

V. RESULTS AND DISCUSSION 

In this section we present our theoretical and simula¬ 
tion results for hard ellipses, freely rotating and orien- 
tationally restricted monolayers of hard ellipsoids. The 
equation of state, the surface coverage, the in-plane and 
out-of-plane order parameters and the isotropic-nematic 
(IN) transition densities are determined for various val¬ 
ues of shape anisotropy. We start with the system of hard 
ellipses, which is obtained by setting the polar angles of 
all ellipsoidal particles {0) to be 7r/2. This condition is 
accompanied by /(cj) = /(0) and that the uniaxial order 
parameter {S) is always -1/2 and the biaxial order param¬ 
eter can be obtained from A = d(j)f{(j)) cos(20). Note 

that A now serves as a 2D orientational order parameter. 
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FIG. 2: Orientational ordering of the system of hard ellipses: 
the equation of state (upper panel) and the two-dimensional 
order parameter vs. density (lower panel). The curves cor¬ 
respond to = 3, 4, 5, 10 and 20 from left to right. The 
curves are the results of PL theory, while the diamond sym¬ 
bols correspond to MC simulation results. The hard ellipse 
system is obtained from oblate ellipsoids by tilting their ori¬ 
entations into the XY plane. 


The theoretical calculations and simulations can also be 
carried out with both oblate-shaped {0 < k < 1) and 
prolate-shaped hard ellipsoids {k > 1), where 0 = nj2 
for all particles. Here we perform the hard ellipse study 
using oblate-shaped ellipsoids. Fig. shows the pressure 
and the 2D order parameter (A) as a function of density 
for = 3, 4, 5, 10 and 20. One can see that more 
particles on the surface (higher densities) are required to 
maintain the same value of the pressure with increasing 
shape anisotropy. This is due to the fact that the XY 
plane becomes more spacious with increasing , i.e. 
a\\ must decrease at fixed a± . 

It can be also seen in Fig. [^that the simulation data 
are well reproduced by the approximate PL theory except 
the values of the order parameters of the N phases cor¬ 
responding to extremely anisotropic particles (/^“^ = 10 
and 20). These data well agree with those given in [33]. 
The order parameter curves reveal for the occurrence of 
isotropic-nematic phase transition, because the phase is 
isotropic at low densities (A = 0), while it is nematic 
at high ones (A > 0). The phase transition is second 
order in the theory, while it is higher order (Kosterlitz- 
Thouless type continuous transition) in the simulation. 



FIG. 3: Shape dependence of the isotropic-nematic (IN) 
transition of hard ellipses: reduced density vs. aspect ra¬ 
tio. The curves are the results of PL theory, while the di¬ 
amond symbols correspond to MG simulation results. The 
inset shows the packing fraction of the IN transition as a 
function of aspect ratio, where r] = p7Ta\\a±/4:. 


Both the theory and the simulation show the same ten¬ 
dencies for IN transition densities and packing fractions 
(see Fig. [^. Making the ellipsoids more anisotropic one 
can see that the IN transition density increases, while 
the IN packing fraction decreases. This apparent con¬ 
tradiction is due to the fact that the increasing shape 
anisotropy makes the system more spacious {a± is con¬ 
stant, while cr|| is decreasing), which requires more par¬ 
ticles for the initiation of the phase transition, while the 
XY plane can be less occupied at the same time. It can 
be also seen that the theory underestimates IN transition 
densities and packing fractions in higher extent with in¬ 
creasing shape anisotropy, which is due to the fact that 
the contribution of higher order virial coefficients are not 
negligible with increasing 

The monolayer of freely rotating hard ellipsoids can 
be achieved by using Oc = 7r/2 and Oc = 0 limiting 
angles for oblate and prolate shapes, respectively. The 
resulting surface coverage, effective aspect ratio, pres¬ 
sure and order parameters are shown for prolate and 
oblate shaped ellipsoids in Figs. |^and[^ respectively. In 
these figures the effective aspect ratio of the correspond¬ 
ing hard ellipse system is given by /^eff = (cr||,e/^±,e) = 

{k/^ sin^ d ^ cos^ 6>), which is the average aspect ra¬ 
tio of the ellipses obtained by intersecting the ellipsoids 
by the XY plane. In turn, the surface coverage can be 
obtained with the help of /^eff as follows: r] = . 

Starting with prolates ellipsoids, one can see that the par¬ 
ticles do not form an isotropic phase, but they align along 
the normal of the XY plane (see the simulation snapshot 
of Fig. even at very low densities. In fact, an isotropic 
phase would appear only at the limit of infinite dilution. 
The density dependence of the surface coverage and the 
effective aspect ratio shows that the particles have less 
and less chance to lie into the XY surface with increas- 
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FIG. 4: Monolayer of freely rotating hard prolate ellipsoids: 
surface coverage (packing fraction) vs. reduced density (up¬ 
per panel) and the equation of state (lower panel) at k = 10. 
The effective aspect ratio of the corresponding hard ellipse 
system (k = (cr\\^e/(^±,e)) is shown in the inset of the upper 
panel, while the out-of-plane order parameter (S) is presented 
in the inset of the lower panel. The curves are the results of 
PL theory, while the diamond symbols correspond to MC sim¬ 
ulation results. The snapshot is the result of replica exchange 
Monte Carlo simulation method. 


ing density because of packing effects. The intersection 
of the ellipsoid with the XY surface is practically a hard 
disk for p* > 0.6, where the effective ellipse aspect ratio is 
almost one. This makes the monolayer of hard ellipsoids 
very similar to the system of hard disks at high densities 
despite the presence of orientational fluctuations. The in¬ 
plane order parameter (Eq. 10) is always zero, while the 
out-of-plane order parameter (Eq. 9) is positive. This 
shows that the out-of-plane orientational ordering is uni¬ 
axial with nematic director parallel to the layer normal. 
It can be also seen that S is very close to its maximum 
value (5'max = 1) at p* =0.7, which corresponds to the 
case where all prolate ellipsoids are parallel and align 
along the layer normal. The reason why prolate ellip¬ 
soids prefer the out-of-plane ordering is that they can 
maximize the free area available on the surface with the 
penalty of orientational entropy loss. Moreover the close 
packing structure of prolate ellipsoids is identical to that 
of 2D hard disks. One can also see that the theory re¬ 
produces quite well the simulation data for all properties 
except the equation of state at high densities. Regarding 
the possible freezing transition at high densities, it has 
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FIG. 5: Monolayer of freely rotating hard oblate ellipsoid 
monolayer: surface coverage (packing fraction) vs. reduced 
density (upper panel) and the equation of state (lower panel) 
at K = 1/10. The effective aspect ratio of the corresponding 
hard ellipse is shown in the inset of the upper panel, while the 
out-of-plane (S) and the in-plane (A) order parameters are 
presented in the inset of the lower panel. The curves are the 
results of PL theory, while the diamond symbols correspond 
to MG simulation results. The snapshot is the result of replica 
exchange Monte Garlo simulation method. 


not been examined here, because it is not the scope of 
our present study. 

The principle of minimizing the intersected area with 
the XY plane applies also for oblate ellipsoids (see Fig. 
1^. The minimal intersected area, which is actually an 
ellipse with cr||^e = ^|| dimensions, can be 

achieved with ordering into the XY plane, while the out- 
plane ordering result in higher intersected area, because 
a\\ < a± for oblate ellipsoids. The resulting in-plane or¬ 
dering can be seen in the simulation snapshot and from 
the high density behavior of the effective aspect ratio 
^ The out of plane order parameter 

goes to -1/2, which corresponds to complete in-plane or¬ 
dering. This in-plane order is isotropic at low densities 
(A = 0), while it is nematic at high densities (A > 0). 
This means that the system undergoes a phase transition 
from a planar uniaxial nematic order (S' < 0, A = 0) to 
a biaxial nematic one (S < 0, A 7^ 0), which can be 
considered as a 2D IN phase transition. The order of 
the uniaxial nematic-biaxial nematic (N-BN) phase tran¬ 
sition is proved to be second order in the approximate PL 
theory, while the transition is higher order continuous in 
the simulation. This transition is the result of the com¬ 
petition between in-plane orientational entropy (favoring 
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FIG. 6: Shape dependence of the uniaxial nematic-biaxial ne¬ 
matic (N-BN) transition of freely rotating hard oblate mono- 
layer: reduced density vs. aspect ratio. The curves are the 
results of PL theory, while the diamond symbols correspond 
to MC simulation results. The inset shows the packing frac¬ 
tion of the N-BN transition as a function of aspect ratio. 


disorder) and the in-plane packing entropy (favoring or¬ 
der). The high density structure of the hard ellipsoid 
monolayer resembles the nematic phase of the 2D hard 
ellipses, while at low densities the structure more or less 
corresponds to an isotropic phase of a polydisperse mix¬ 
ture of hard ellipses. One can also see that the agreement 
between the theory and the simulation is quite good for 
oblate ellipsoids, especially with respect to the equation 
of state. An interesting feature is that the theory overes¬ 
timates the out-of-plane ordering; this could stem from 
the effect of higher virial coefficients, which are not in¬ 
cluded in the theory. 

Fig. 1^ shows the in-plane isotropic-nematic (or N-BN) 
transition densities as a function of aspect ratio result¬ 
ing from the PL theory and simulation. On the other 
hand, the transition density increases with increasing 
shape anisotropy, while the packing fraction decreases. 
This is due to the fact that we go to the "volumeless” 
hard needle limit with decreasing k , where the transi¬ 
tion density saturates at finite value, while the packing 
fraction vanishes (note that /^eff ^ 0 and 77 = p* 7 r/^eff/ 4 ). 
The agreement for the transition densities is again very 
good between the theory and simulation. To see the effect 
of out-of-plane orientational freedom it is worth plotting 
the IN densities of 2D hard ellipse system and those of 
freely rotating hard oblate monolayer together (see Fig. 
[^. It can bee seen that the simulation does not show 
a substantial difference between the two systems, while 
the theory predicts that the out-of-plane freedom desta¬ 
bilizes the nematic order, i.e. the IN density curve is 
shifted into the direction of higher densities. In the sim¬ 
ulation the nematic phase evolves from a very ordered 
planar phase (isotropic in the XY plane), where the ef¬ 
fective aspect ratio is almost identical with the aspect 
ratio of the ellipsoids (see Fig. [^, i.e. the correspond¬ 
ing hard ellipse system is almost monodisperse. This is 


FIG. 7: Gomparison of the transition densities of freely ro¬ 
tating hard ellipses and confined hard oblates in the reduced 
density-aspect ratio plane. The continuous (ellipsoids) and 
dashed (ellipses) curves are the results of PL theory, while 
the open and filled diamond symbols correspond to MG sim¬ 
ulation results for hard ellipses and hard ellipsoids, respec¬ 
tively. The inset shows the packing fractions of IN and N-BN 
transitions as a function of the aspect ratio. 


not the case in the theory, where the transition happens 
at such densities, where the corresponding hard ellipse 
system is still polydisperse with larger shape anisotropy 
(/^eff) than . As a result the effect of out-of-plane orien¬ 
tational freedom is more pronounced in the theory than 
in the simulation. This shows that the theory exagger¬ 
ates the effect of out-of-plane orientational freedom on 
the N-BN phase transition. 

Finally we show how the IN transition properties 
change with switching on the out of plane orientational 
freedom through the varying limiting polar angle (Oc)^ 
which is between 0 and 7 r /2 for both prolate and oblate 
shapes. If Oc = 7 r /2 (cos 6 >c = 0) the prolate ellipsoids 
on the XY plane behave like the system of hard ellipses, 
i.e. they undergo a 2D IN phase transition with increas¬ 
ing density (see Fig. |^. However this phase transition 
is destabilized with decreasing the limiting angle (Oc)^ 
which corresponds to increasing cos Oc- As the orienta¬ 
tional window is widened (decreasing 6 >c), the prolate el¬ 
lipsoids tend to minimize their occupied area on the XY 
surface, which results in less anisotropic in-plane ellipses 
and higher transition densities for all studied aspect ra¬ 
tios. In addition to this, the in-plane order transforms 
continuously into out-of-plane order, i.e. S becomes posi¬ 
tive with cos Oc (see Fig. [^. Hence, the 2D IN transitions 
turn into N-BN transitions for confined prelates which 
are lost for the free rotating case. The surface cover¬ 
age (packing fraction) curves show very clearly that the 
occupied area on the XY surface increases enormously 
with decreasing the polar angle restriction, i.e. the 2D 
IN phase transition is destabilized with decreasing Oc . 
The value of Oc where the orientational ordering transi¬ 
tion is preempted by the positional one (freezing transi¬ 
tion) cannot be determined with the present PL theory. 









Interestingly, MC simulation shows that the IN density 
is practically not affected by the value of limiting polar 
angle, but the transition disappears at values of Oc close 
to zero. Note that the IN transition density curves move 
into the direction of lower densities with increasing shape 
anisotropy because the occupied area of the ellipsoids be¬ 
comes larger with increasing k as > aj_. The case of 
oblate ellipsoids is different because Oc = 7r/2 corresponds 
to the freely rotating limit, while = 0 is the hard disk 
limit. Starting from the freely rotating case one can see 
that the IN transition density shows a maximum at an 
intermediate value of Oc , which can be attributed to the 
combined effect of decreasing in-plane shape anisotropy 
and the increasing intersected area of the ellipsoids with 
decreasing Oc • The N-BN transition packing fraction be¬ 
haves simply since it is always an increasing function of 
cos 6>c, i.e the nematic phase is destabilized. This means 
that the nematic-biaxial nematic curve and the isotropic- 
solid curves must cross each other at a threshold value of 
cos Oc making the nematic order metastable and produc¬ 
ing a plastic solid (an orientationally disordered solid). 
This plastic solid would eventually turn into an ordered 
solid at higher densities. However, this value of cos Oc has 
not been searched by the simulation and theory, because 
our study focuses only on the orientational ordering. Fi¬ 
nally, the S order parameter shows that the system un¬ 
dergoes a structural change from in-plane order into out- 
of-plane order at the N-BN transition. Our results show 
that the effect of out-of-plane orientational freedom is to 
minimize the occupied area on the surface which allows 
the most efficient packing of the ellipsoids on the XY 
plane. 


VI. CONCLUSIONS 

We have examined the orientational ordering proper¬ 
ties of the monolayer of hard ellipsoids using replica ex¬ 
change Monte Carlo simulation method and the Parsons- 
Lee density functional theory. We have found that both 
the shape anisotropy and the orientational restriction af¬ 
fect substantially the orientational order on the plane. 
The subtle interplay between orientational and packing 
entropies results in different structures and phase be¬ 
haviors for prolate and oblate shaped ellipsoids. The 
main driving force of the ordering behavior is to mini¬ 
mize the intersected area between the particle and the 
surface to realize the close packing structure with in¬ 
creasing surface density. The minimal intersected area 
is a disk for prolate, while it is an ellipse for oblate el¬ 
lipsoids. Highly packed structures can be achieved with 
only out-of-plane (in-plane) ordering for prolate (oblate) 
ellipsoids, i.e. prolate ellipsoids prefer to order along the 
normal of the plane, while oblate ellipsoids like to order in 
the plane. The monolayer of prolate ellipsoids with very 
small out-of-plane orientational freedom behaves almost 
identically to the system of 2D hard ellipses, i.e. it form 
2D isotropic and nematic phases. The gradual rise of the 



FIG. 8: Monolayers of orientationally restricted hard ellip¬ 
soids: N-BN transition density vs. cos^c- The upper panel 
is for prolate shapes while the lower one for oblate shapes. 
Insets show the cos Oc dependence of the surface coverage and 
the out-of-plane order parameter at the N-BN ordering tran¬ 
sition. The curves are the results of PL theory. 


out-of-plane freedom, through decreasing Oc , allows the 
ellipsoid particle to lean out from the plane and to de¬ 
crease the intersected area with the plane. This involves 
less anisotropic shape in the interactions and destabi¬ 
lization of the biaxial nematic phase (2D nematic) with 
respect to uniaxial nematic one (2D isotropic). At full 
orientational freedom {Q < 0 < it) prolate ellipsoids do 
not form an isotropic phase, but they are ordered along 
the normal of the confining plane even at very low densi¬ 
ties and behave similarly to the 2D system of hard disks 
at high densities. 

The monolayer of oblate ellipsoids without out-of-plane 
freedom is identical to the 2D system of hard disks. The 
gradually increasing freedom in polar angle {0) allows 
the oblate ellipsoid to decrease its intersected area with 
the confining plane through leaning out from the plane. 
This makes the ordering planar for the main symmetry 
axis of the ellipsoid. However the anisotropic interactions 
between the tilted ellipsoids give rise to an additional in¬ 
plane order, i.e. the phase is biaxial nematic, which is 
due to the excluded area gain coming from the in-plane 
ordering of elliptical intersections. In the freely rotating 
case the oblate ellipsoids form a planar nematic phase 
at low densities and a biaxial nematic one at high densi¬ 
ties. Interestingly, the IN density of 2D ellipses is almost 
identical to that of freely rotating ellipsoids confined to 













9 


a plane. In summary the out-of-plane orientational free¬ 
dom stabilizes the in-plane nematic ordering in the mono- 
layer of oblate ellipsoids, while the opposite happens in 
the monolayer of prolate ellipsoids. This is congruent 
with the 2D ellipse and 2D disk phase diagrams, which 
are obtained as limiting cases in this study. In our sim¬ 
ple model we have neglected the effect of out-of-plane 
positional freedom which we expect to have a relatively 
small effect on the system behavior. On the other hand, 
not including a soft wall-particle interaction may lead 
to deviations from experimental set-ups where capillary 
phenomena are frequently important. No doubt, the in¬ 
clusion of this interaction would substantially increase 


the computation burden of the problem. We leave this 
issue for future studies. 
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